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Abstract: The coupling effects between the mechanical and electric properties of 
piezoelectric materials have drawn significant attention for their potential applications as 
sensors and actuators. In this investigation, two piezoelectric actuators are symmetrically 
surface bonded on a cross-ply composite laminate. Electric voltages with the same 
amplitude and opposite sign are applied to the two symmetric piezoelectric actuators, 
resulting in the bending effect on the laminated plate. The bending moment is derived by 
using the classical laminate theory and piezoelectricity. The analytical solution of the 
flexural displacement of the simply supported composite plate subjected to the bending 
moment is solved by using the plate theory. The analytical solution is compared with the 
finite element solution to show the validation of present approach. The effects of the size 
and location of the piezoelectric actuators on the response of the composite laminate are 
presented through a parametric study. A simple model incorporating the classical laminate 
theory and plate theory is presented to predict the deformed shape of the simply supported 
laminate plate. 
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1. Introduction 

Piezoelectric materials with the advantages of quick response, low power consumption and high 
linearity have drawn much attention in the past decade. Piezoelectric devices are of great interest in 
structural engineering with applications to shape control, vibration suppression and noise reduction [1]. 
Smart structures integrated with sensors and actuators have the capability to respond to a changing 
environment and control the structural movement. Piezoceramics are the most common material used 
in smart structures and can be surface bonded to existing structures to form an online monitoring 
system, or embedded in composite structures without significantly changing the structural stiffness 
system. Bailey and Hubbard [2] developed first adaptive structure using polyvinylidene fluoride 
(PVDF) film as actuators to control the structural vibration of a cantilever beam. Crawley and 
de Luis [3] studied a beam with surface bonded and embedded piezoelectric actuators to investigate 
the load transfer between the actuator and host beam. Huang and Sun [4] studied the load transfer and 
wave propagation of an anisotropic elastic medium induced by the surface bonded piezoelectric 
actuator. Dimitriadis et al. [5] used two dimensional patches of piezoelectric material bonded to the 
surface of a simply supported plate as vibration actuators to excite the selected modes. Shape control is 
one of the major applications for piezoelectric materials. Koconis et al. [6] controlled the shape of 
composite plates and shells with piezoelectric actuators. Luo and Tong [7] developed a finite element 
model to simulate twisting and bending shape control using the ortho tropic piezoelectric actuators. Lin 
and Nien [8] used piezoelectric actuators to control the deflection and shape of the composite 
laminates. Other works related to the shape control of composites with piezoelectric actuators were 
presented by Bulter et al. [9,10]. The finite element method, a widely accepted and powerful tool for 
analyzing complex structures, is capable of dealing with the piezoelectric smart structures. Numerous 
studies have been completed on analyzing piezoelectric structures to account for the piezoelectric 
effect [11-13]. Actuators used in active control of smart structures have to be located appropriately to 
ensure maximum control and measurement effectiveness. The methods for optimal placement of 
sensors and actuators were investigated by many researchers [14-16]. 

The present work investigated the load transfer between surface bonded piezoelectric actuators and 
the host structure. The proposed method is an extension of the one dimensional beam embedded with 
PZT derived by Crawley and de Luis [3] and two dimensional plate surface bonded with PZT derived 
by Dimitriadis et al. [5]. The model consists of two piezoelectric actuators symmetrically surface 
bonded on a cross-ply composite laminate subjected to electrical voltage. An analytical expression of 
the bending moment induced by the piezoelectric actuators was derived by incorporating the 
composite laminate theory and piezoelectric effect. The bending moment was then applied to the 
composite laminate with simply supported boundary conditions. A closed form solution to the 
deflection of the simply supported composite laminate was obtained by using the plate theory. The 
analytical solution was validated with the finite element results. Finite element method has been 
widely used in structural analysis. It is reasonable to verify present approach with finite element 
method. The effects of the size and location of the piezoelectric actuators on the response of the 
composite laminate are presented through a parametric study. The objective of this investigation is to 
develop an analytical expression of the response of a thin plate excited by the bonded piezoelectric 
actuators. A simple model incorporating the classical laminate theory and plate theory is proposed to 
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predict the deformation of the composite laminate plate. The feasibility of controlling the deflected 
shape of the plate is illustrated by placing the actuators at various locations. 



2. Bending Moment 



Consider two piezoelectric actuators symmetrically surface bonded on a cross-ply composite 
laminate. The polarized direction is along the z-axis. For an unconstrained thin piezoelectric actuator, 
equal strains in both x and y directions will be induced when activated by a voltage along the poling 
direction. The magnitude of the strain can be expressed in terms of the piezoelectric constant dn, 
applied voltage Vand actuator thickness t pe , as follows: 

{s x ) pe ={s y ) pe =s pe =^V (1) 

V 

where subscripts pe and p represent the quantities associated with piezoelectric actuator and host plate, 
respectively, throughout this paper. When an electrical field is applied in the direction normal to the 
actuator surface, surface strains are generated Equation (1). Due to the coupling of the actuator to the 
structure, forces and moments are induced in the bonded area of the structure. Since this work focus on 
the deformation of the plate induced by the bending moment, only constant d31 is considered in 
this model. 

The two actuators are activated by applying a voltage of equal magnitude and opposite sign to the 
opposing actuators. The opposite directions of the surface tractions at the interfaces between the 
actuator and plate cause the uniform bending moments along the actuator boundaries as shown in 
Figure 1. 

Figure 1. Bending moment acting on the composite laminate induced by the lead zirconate 
titanate PZT actuators. 
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The symmetry of the actuators with respect to the midplane (z = 0) results in no net extension or 
contraction in the midplane of the plate. In the following derivation, the piezoelectric actuators are 
assumed to be perfectly bonded to the composite laminate and in the state of plane stress. This implies 
the strain continuity across the interfaces as shown in Figure 2. 

Figure 2. Strain distribution across the thickness of the composite laminate. 



Piezoelectric 




The strains across the thickness of the composite laminate due to the bending moment can be 
expressed as follows: 

s x = Zk x ; s y = Zk y ; s xy = Z/c xy (2) 

where Z is the position along the thickness measured from the midplane of the laminate, k x , K y and 
K xy are the curvatures of the composites laminate. 

The bending stresses in the k-th layer of the composite laminate are: 
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where E\ and Ei are the Young's moduli along the fiber direction and normal to the fiber direction, 
respectively. Gn is the shear modulus, v i2 and v 2 i are the Poisson's ratios, 9 is the fiber orientation 
measured from the global x-axis. 

In the case of cross-ply (9 = 0° or 90°), Q l6 and Q 26 are equal to zero. The bending stresses in the 

k-th layer can be reduced to: 

<t™=Z&Pk x+ q£W (4a) 

^=z[Q^K x+ Q^ Ky \ (4b) 

< ] = ZQ^K xy (4c) 

The bending stresses in the piezoelectric actuator are: 

E 



pe 1 - v 

pe 



(5a) 
(5b) 



pe 



where E pe and v pe are the Young's modulus and Poisson's ratio of the piezoelectric actuator. 

The bending moments per unit length m x and m y are defined as the forces Oxdz and Oydz times the 
moment arm z, respectively. The sum of the bending moments with respect to the neutral axis 
(z = 0) across the thicknesses of the composite laminate and piezoelectric actuators are zero, such that: 

f (<r x ) ZdZ + {a x ) ZdZ + JT* {a x ) ZdZ = 0 (6a) 

•>— r p J—n—t P e *t P e 

i (°> ) p zdz + el k ^ + n t zjz = ° (6b) 

where t is one half of the composite laminate thickness and h is the thickness of the 
piezoelectric actuator. 

l(Ai), +2(Ai) pe K + l D n) p +2(A 2 )^K ={ l + v peH B n)peZpe (7a) 

l(A 2 )„ + 2(A 2 R + [(A 2 ), + 2(D 22 K=( 1 + v pe M*i i U ^ ( 7b ) 
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where and Z k .\ represent the positions of the top and bottom surfaces of the k-th layer in the 

composite laminate, respectively. 

The curvatures k x and K y can be solved from Equation (7) as: 

K x =A£pe (8a) 

K y= A 2^ P e (8b) 

A = B npe (\ + v pe ) | (((A 2 ) P +A2pJ-^ e ((Ai)p+AipJl + v p J+^np^(A2)p+A2pJl + Vp e ))) (8c) 
1 (D n ) p+ D Upe (((D u ) p +D Upe \-(( D n) P +Dn pe f+(( D n) P + D npe \{D 22 ) p + D 22pe ))) 

a = -B Upe ((D u ) p + D Up J\ + v pe )+B npe ((D l2 ) p + D npe \\ + v pe ) 

2 "((A 2 )p + D nJ + +D Up J(D 22 ) p+ D 22pe ) 

The bending moments per unit length m x and m y acting on the composite laminate as shown in 

Figure 2 can be calculated: 

^ = t€ Z(QS^+QS\W (9a) 
=£[C ^>] = t€ Z(Q{?K x+ Q^ y )d Z ] (9b) 

fc=l *=1 

Substituting Equation (8) into Equation (9), leads to the bending moments: 

m x =C x s pe (10a) 

m y =C 2 s pe (10b) 

c^a^d^+a^d^) 

C 2 =(A 1 (D 12 ) p +A 2 (D 22 ) p ) 



3. Deflection of a Simply Supported Composite Plate 



A rectangular composite laminate plate with simply supported boundary conditions is considered in 
this work. The location of the surface bonded actuator is shown in Figure 3. The activated piezoelectric 
actuators will induce bending moments as derived in Equation (10) to the composite plate and can be 
expressed in terms of unit step functions as follows: 

m x =C l £ pe [h(x-x l )-h(x-x 2 )][h(y-y l )-h(y-y 2 )] (11a) 
m y = C 2 £ pe [h(x-x 1 )-h(x-x 2 )][h(y-y 1 )-h(y-y 2 )] (lib) 
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Figure 3. Surface bonded actuator on the composite laminate plate. 




•Piezoelectric 



Using the classical plate theory, the equilibrium equation for the plate can be written in terms of the 
plate internal moments M x ,M y ,M xy and the actuators induced moments m x ,m y as: 



d\M x -m x ) , 2 8 2 M xy | d 2 (M y -m y ) =Q 



(12) 



dx z dxdy dy 1 

The internal moments M x ,M y ,M xy can be expressed in terms of the flexural displacement w . 
Moving the moments m x ,m y to the right hand side of Equation (12), yields to the following 
equilibrium equation of the composite laminate plate: 
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(13) 



P = 



■ + ■ 



dx 2 dy 2 

where (D ll ) p ,(D 22 ) p ,(D 66 ) p are the bending stiffness of the composite laminate as defined in 
Equation (7). 

Substituting Equation (11) into Equation (13), leads to the governing differential equation: 



d 4 w 



d 4 w 



d 4 w 
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(14) 



For a simply supported rectangular plate, the flexural displacement w can be expressed by the 
following Fourier series: 



M.x,y) = YJTw m 



. mwc . nm 
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Substituting Equation (15a) into Equation (14), solve for the constant W mn as follows [5]: 



W = 

mn 4 4 



2 2 2 2 4 4 

mn' . mn nn - TT n n 
— r(Ai)» + — o 7^ 2H x +—^(^22), 



(15b) 
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l cos r,„*i - cos r m x 2 k cos r„yi - cos 
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a b 



(15d) 



4. Finite Element Analysis 



The finite element method is a widely used and powerful tool for analyzing complex structures. It is 
capable of dealing with the piezoelectrical materials. Many researchers have modelled the 
piezoelectric actuation using the finite element method. A commercially available finite element 
software ANSYS has the ability to analyze the piezoelectrical materials. In this study, the finite 
element software ANSYS is adopted to investigate the deflection of a simply supported plate induced 
by the surface bonded piezoelectric actuators. To perform the ANSYS finite element analysis for the 
piezoelectric actuator bonded structure, SOLID 45 elements and SOLID 5 were used in the composite 
plate and piezoelectric actuators, respectively. 

A typical three dimensional finite element mesh is shown in Figure 4. A voltage between the upper 
and lower surfaces of the SOLID 5 elements is applied, which results in an electric field along the 
poling direction of the actuator. The deflections obtained from the finite element method are compared 
with the analytical solutions of Equation (15) to validate the present approach. 



Figure 4. 3-D finite element mesh. 
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5. Numerical Validation and Examples 



In the following numerical examples, the composite material is carbon/epoxy with stacking 
sequence [0/90/90/0]. The composite material properties of carbon/epoxy are listed in Table 1. The 
dimensions of the composite laminate plate are length a = 0.38 m, width b = 0.3 m, thickness 
t p = 1.5876 mm. The piezoelectric actuator is assumed to be PZT G-1195 with the material 
properties [17] of Young's modulus E pe = 63 GPa, Poisson's ratio Vp e =0.3, density p pe = 7,600 kg/m , 
piezoelectric constant dn = 1.9 x 10~ 10 V/m and thickness t pe = 0.15876 mm. The effects of the size 
and location of the actuators are presented through a parametric study to investigate the deflection and 
deformed shape of the composite plate activated by the surface bonded piezoelectric actuators. 



Table 1. Material properties of carbon/epoxy. 



Longitudinal 


Transverse 


Shear 


Shear 


Poisson's 


Poisson's 


modulus Ei 


modulus E2 


modulus G12 


modulus G23 


ratio V12 


ratio V23 


108 GPa 


10.3 GPa 


7.13 GPa 


4.02 GPa 


0.28 


0.28 



5.1. Example 1: Three different sizes of actuators 



Two piezoelectric actuators are surface bonded on the top and bottom surfaces of the composite 
plate. Three different sizes of piezoelectric actuators with the dimensions of 0.06 m x 0.04 m, 
0.08 m x 0.06 m and 0.1 m x 0.08 m, respectively, bonded on the central area of the composite plate as 
shown in Figure 5 are considered in this example. The voltages of +1 V and - IV are applied to the top 
and bottom actuators, respectively, resulting in a bending moment acting on the composite plate. The 
deflections of the composite plate induced by the actuators are calculated using both the analytical 
prediction of Equation (15) and finite element method. 

Figure 5. Three different sizes of PZT actuators. 
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Figure 6 shows the deformed shape of the composite plate predicted by Equation (15) and finite 
element method. The flexural displacements along the horizontal line y = b/2 and the vertical line 
x = all of the composite plate are presented in Figure 7. The deflection is increasing as the size of 
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actuators increases. The deflections of the composite plate predicted by Equation (15) and finite 
element method are agree well. The maximum deflections in the composite plate induced by the three 
different sizes of actuators are listed in Table 2. It shows that the difference between the present 
approach and finite element method is within 7 %. Similar difference was observed by Qing et al. [18] 
as comparing their semi-analytical solutions with ANSYS finite element results. The finite element 
methods can handle numerically a lot of problems in engineering field and has proved to be a powerful 
tool for the design and analysis of piezoelectric devices. Owing to the boundaries and geometric 
complexity of the adaptive structure, it is reasonable to validate the present model with finite 
element method. 

Figure 6. Deformation of the composite plate induced by the PZT actuator 
0.06 m x 0.04 m obtained by (a) Equation (15) (b) finite element method (unit: m). 
(a) Deformed shape obtained by Equation (15). (b) Deformed shape obtained by ANSYS. 

(a) (b) 




Figure 7. Displacements of the composite plate induced by PZT actuators with three 
different sizes. 
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Figure 7. Cont. 
(b) Along the vertical line x = all. 




□ 


0.06*0.04m ANSYS 




— 0.06*0.04m eq(15) 




0. 08 *0.06m ANSYS 




— 0.08*0. 06m eq(15) 


o 


0.1*0. 08m ANSYS 




- 0.1*0. 8m eq(15) 



(b) vertical line x=a/2 (m) 



Table 2. maximum deflection of the composite plate induced by the PZT actuators with 
three different sizes. 



^^^^lethod 
Size ^^^^ 


ANSYS 


Equation (15) 


Error 

(%) 


PZT 
0.06 m x 0.04 m 


1.68xl0" 3 mm 


1.59xl0~ 3 mm 


5.4 


PZT 
0.08 mx 0.06 m 


2.90xl(T 3 mm 


2.73xl0~ 3 mm 


6.3 


PZT 
0.1 mx 0.08m 


4.26xl0" 3 mm 


3.98xl0~ 3 mm 


6.8 



5.2. Example 2: Three different locations of actuators 

To study the capability of control the deflection shape of the plate, actuators are placed at various 
locations. In this example, the piezoelectric actuators are surface bonded at three different locations, 
central, right and top region of the plate, respectively, as shown in Figure 8. These three typical 
locations were arbitrarily chosen to demonstrate the influence of the actuator location on the deflection. 
The deformed shapes of the composite plate induced by the PZT actuators surface bonded at the 
central, right and top region of the plate are shown in Figures 6, 9 and 10, respectively. The difference 
of the deflected curves shown in Figures 6, 9 and 10 demonstrates that the shape of the plate can be 
controlled by placing the actuators at various locations. The flexural displacements of the plate along 
the horizontal line y = b/2 and vertical line x = a/2 are presented in Figure 1 1 . The deflections of the 
plate obtained by the present approach of Equation (15) and finite element method are in close 
agreement. Table 3 lists the maximum deflections of the plate induced by the piezoelectric actuators 
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surface bonding at three different locations. It shows that the difference between the present approach 
and finite element method is within 8%. 

Table 3. maximum deflection of the composite plate induced by the PZT actuators at three 
different locations. 



Location^\ 


ANSYS 


Equation (15) 


Error 

(%) 


PZT at 
central region 


1.68xl0" 3 mm 


1.59x10 3 mm 


5.4 


PZT at 
right region 


8.62x10^ mm 


8.00xl(T 4 mm 


7.7 


PZT at 
top region 


1.27xlCT 3 mm 


1.34xl0" 3 mm 


4.9 



Figure 8. Three different locations of the actuator. 
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Figure 9. deformation of the composite plate induced by the PZT actuator bonded on the 
right region calculated by (a) Equation (15) (b) finite element method, (a) Deformed shape 
calculated by Equation (15). (b) Deformed shape calculated by ANSYS. 
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Figure 10. deformation of the composite plate induced by the PZT actuator bonded on the 
top region calculated by (a) Equation (15) (b) finite element method, (a) deformed shape 
calculated by Equation (15). (b) deformed shape calculated by ANSYS. 




Figure 11. Flexural displacement of the composite plate obtained by ANSYS and 
Equation (15) along the horizontal line (y = b/2) and vertical line (x = a/2) for three 
different locations of actuator. 
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6. Conclusions 

Piezoelectric materials are often used as strain actuators and shape control of smart structures, as 
they are compact and response quickly. In this investigation, two piezoelectric actuators are 
symmetrically surface bonded on a composite laminate plate. Electric voltages with the same 
amplitude and opposite sign are applied to the two symmetric piezoelectric actuators, resulting in the 
bending effect on the plate. Theoretical model of the bending moment is derived by using the theory of 
elasticity to represent the interaction of the actuator and the host plate. Following the classical plate 
theory, the deflection of a simply supported plate subjected to the bending moment can be obtained. 
The analytical solutions are validated with the finite element results. The effects of size and location of 
actuators on the responses of the plate are presented through the parametric study. Utilization of the 
laminate theory and plate theory, the deformed shape of the laminate plate can be predicted 
analytically. The methodology proposed in this paper is easy to employ, and provides an alternate way 
of solving this complicated problem analytically with accurate results. 
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